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Abstract. Neutrino oscillations have become well-known phenomenon; the 
measurements of neutrino mixing angles and mass squared differences are continuously 
improving. Future oscillation experiments will eventually determine the remaining 
unknown neutrino parameters, namely, the mass ordering, normal or inverted, and the 
CP-violating phase. On the other hand, the absolute mass scale of neutrinos could 
be probed by cosmological observations, single beta decay as well as by neutrinoless 
double beta decay experiments. Furthermore, the last one may shed light on the 
nature of neutrinos, Dirac or Majorana, by measuring the effective Majorana mass of 
neutrinos. However, the neutrino mass generation mechanism remains unknown. A 
well-motivated phenomenological approach to search for new physics, in the neutrino 
sector, is that of non-standard interactions. In this short review, the current constraints 
in this picture, as well as the perspectives from future experiments, are discussed. 
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1. Introduction 

The last two decades have seen the great success of many neutrino experiments, in 
particular, those that have contributed to the first observations of the various types 
of neutrino oscillation m 121 El m 13 El El- The “standard” oscillation parameters, 
defined in the framework of the three active neutrinos, have been determined with 
significant accuracy [8] apart from the yet unknown mass ordering, or the sign of the 
mass squared difference relevant for oscillation of atmospheric neutrinos, as well as 
the CP-violating phase, the most elusive parameter. Non-oscillation experiments that 
have measured the neutrino cross section with high accuracy, also provided valuable 
information for neutrino physics lama ITT] . On the other hand, there are oscillation 
experiments H21 EH 0 UM that have observed hints for neutrino oscillation into an 
additional sterile neutrino state. 

So far, there is no experimental evidence that neutrinos pose some non-standard 
properties beyond masses and mixing or some extra new interactions, different from the 
weak interaction, not described by the Standard Model (SM). Such interactions, often 
called non-standard interactions (NSI) of neutrinos, if they exist, are interesting from a 
phenomenological point of view, since they directly indicate the presence of some new 
physics beyond SM. 

Possible presence of NSI was first pointed out by Wolfcnstein [[16j □n, followed 
by the works done in [18, T9j 120! [2T, 22], [23] [21, [25] in the early stage (before the 
experimental discovery of neutrino oscillation m) and afterwards, studied by a large 
number of authors, for example, in [26] [271 ESI ESI ESI [3TJ [32], [33] [3H [35] [36] [371 EH ESI 

SO] SI] S21 S3] SH SHI SHI HZ] SHI SSI E01 EH E21E3] EH ESI EH] EZ] [3H1 ESI EH] ED E21 ESI 

EH [65] [66] [67] ESI [69] [70] [TI] [72] [73] EH [75] [76] [TT] [78] [79]. In this short review, the 
current status of NSI is discussed, mainly from the phenomenological point of view. 

Phenomenologically, NSI can be described with an effective four fermion 
Lagrangian SSI, 

- C‘Ji = ei P e 2V2G F (f alp L^)(fYPf ), (1) 

where Gp is the Fermi constant, is the parameter which describes the strength of 
the NSI, / is a first generation SM fermion (e,w or d), P denotes the chiral projector 
(L, R = (l±7 5 )/2), and a and fd denote the neutrino flavors: e, /i or r. This Lagrangian 
describes neutral current (NC) interactions and they will be the focus of this work. 

While the Lagrangian in ([!]) provides the general description of NSI, it is also 
possible to study NSI using other approaches, defining the parameters depending on 
whether we are considering the neutrino at the production point, e , during propagation 
(taking into account matter effects) e m , or at the detection point, e D . The reader must 
be aware that there is no universal notation for the NSI parameters, and some authors 
use similar notation for different quantities, although the definition of the parameters 
in terms of Eqs. ([!]) and (j2j) (see below) is more standard. 

Let us consider the presence of NSI in the neutrino source with more detail. In 
general, source neutrino fluxes are produced from charged current (CC) interactions such 
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as pion and muon decay and, in the presence of additional non-standard CC interactions, 
it would be necessary to consider contributions that would include terms proportional 
to 


2G f - 75 ) 7 ^]- 

a 


( 2 ) 


In this case, the experimental value of the Fermi constant will be given by T(i] 
G exp = Gp\J |1 + £^ e c \ 2 + \£° c \ 2 + kT C l 2 an( i wou ld imply that the new interactions 
could be parameterized as [3Tj 
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A similar expression could be obtained for the case of neutrino detection e D if we 
are interested in CC interactions. It is important to notice that for NC NSI the 
corresponding expression for e s ’ D is different. As will be discussed below, for NSI coming 
from NC, the left and right couplings e L,R appears naturally, while the e^m for the CC 
case are considered to be left handed. 

For the case of neutrino propagation, there is a direct relation between e m and the 
NSI parameters coming from the Lagrangian in Eq. (JTj) . It will be seen in section [ 3 ] 
that, during propagation in matter, the neutrino potential will be sensitive only to 
vector currents (C = e L + e R ) and, therefore, we will end up with the relation, 
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where Vf = V/(r) = \Z2G F Nf(r) with / = e,w or d (see section 3). 

Solar, atmospheric and long baseline neutrino oscillation experiments are expected 
to give better constraints on propagation NSI parameters coming from matter effects, 
while non oscillation experiments are more sensitive to NSI in production and/or 
detection. Both types of experiments provide valuable complementary information on 
NSI. One of the main disadvantages of non-oscillation experiments is that the flavor 
changing NSI will be present in the interaction only at the second order level (£ 2 ) while 
propagation effects appear at first order (e). However, non-oscillation experiments also 
have some advantage, for instance, they may also be sensitive to axial currents, while 
oscillation experiments are not. 

A model independent analysis that considers all the contributions coming from 
Eq. (jlj) will imply a large number of free parameters, e{^. To our knowledge, an analysis 
considering all the NSI contributions at the same time has never been done. I 11 practice, 
one must constrain the analysis to a number of parameters that could be handled by 
current computational methods and give useful information about the freedom for new 
physics in the neutrino sector. 

In this work we will discuss the current status of NSI studies. We will start by 
giving some examples, in section [2j of the kind of new physics that can be tested by 
using this formalism. Although we will stress the case of NC interaction, the result 
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could be converted into constraints for new physics in the CC sector. We will show 
the current constraints coming from both oscillation and non-oscillation experiments 
in sections [3] and [4] along with a brief explanation of the phenomenological procedure 
to obtain such bounds. Future perspectives to improve the current constraints will be 
discussed in section [5j Finally, conclusions will be given in section |6j 

2. NSI and models for new physics 

Although the NSI formalism appears as a correction to the vector and axial couplings, 
it can account for different types of new physics. In this section, three different classes 
of Standard Model extensions are described, in term of the NSI parameters, and a 
particular example is shown in every case as an illustration of the formalism. 


2.1. Extended gauge symmetries 


Any extension of the SM local gauge symmetry SU(2)l ® U (l)y, in general, introduces 
new gauge bosons that modify the vector and axial coupling. Typical examples are E§ 
string inspired models, that at low energies introduce the extra groups U{l) x ® U{ 1)^, 
leading to an additional Z' neutral gauge bosons. A similar situation happens with the 
left-right symmetric model SU {2) L ® U(1)y ® SU(2) R ®U(l)'y where one extra neutral 
and one extra charged gauge bosons, Z' and IF', appear. 

Within the NSI formalism, there is an easy direct relation between the 
phenomenological parameters and the parameters coming from these models. We can 
consider, for instance, the neutrino dispersion off nuclei, described in the SM Lagrangian 

Ci? = -% E MU - 7>J {/ aL [?7„(1 - 7 5 )«] + f R Ml + 7 5 )9] } . (5) 

v q=u,d 

where f qL ’ R are the SM coupling constants defined elsewhere [8j. For the case of E 6 
models, where two additional neutral vector bosons arise, there will be an additional 
contribution to these couplings constants given by [TS] 


^ = _ 4 M 

Ml 
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( 6 ) 


where Mz is the mass of the SM neutral gauge boson; Mz> accounts for the mass of an 
additional, heavier, new gauge boson; 9w is the weak mixing angle; p^ is the parameter 
which accounts for the radiative corrections; and the angle f3 describes the mixing 
between the two extra gauge bosons that arise from the U( l) x and U(l)^ symmetries. 
These models have yet another extra gauge boson that is considered to be heavier than 
the Z' and decoupled from the above Lagrangian. 










NSI 


5 


2.2. Additional neutral leptons 

An extension of the SM fermion content can give rise to a rich phenomenology, especially 
when it contains extra neutral leptons, usually assumed to be heavy. Within this 
framework, it is possible to work in the standard SU( 2) (8) 1/(1) gauge symmetry and 
consider the additional mixing of isodoublet and isosinglet neutral leptons [80]. As a 
result, the V — A couplings will deviate from the SM prediction. For example, for the 
case of the CC, the ordinary light isodoublet neutrinos will mix with the extra heavy 
isosinglets; this mixing will be described by a matrix 

K=(K l ,K h ), (7) 

where K L and K H describe, respectively, the mixing of ordinary isodoublet light 
neutrinos and extra heavy isosinglets. Notice that the signals of new physics coming from 
this matrix might appear in charged currents through the deviations from the standard 
interactions (that is, in the detection) as well as in oscillation experiments, since the 
mixing matrix Kj j is no longer unitary. Moreover, regarding the NC interactions, they 
will be described by the interaction 

C = 1 r^ ?r Z, v - Lll ,K^K„ L , ( 8 ) 

Z Sill o\y 

where the matrix K^K can be considered as a natural source for NSI in the neutral 
sector [80] . 

The enriched structure that arise from these models can be parameterized in 
different forms that must take into account all the new mixing angles and phases. One 
of the most studied schemes in this context is the seesaw model [STJ B21 [H3, ST) [SO] , 
which gives a natural explanation for the smallness of the neutrino mass. In these 
models, however, the sizeable signals at low energies are expected to be negligible and, 
therefore, the effective NSI should be negligible. There are, however, other models where 
the low energy effects, though small, may be sizeable in the near future, such as the 
so-called inverse seesaw model [SSI E 5 j • 

Although it is not common in the literature to consider these analyses in term of 
the NSI formalism, it is possible to include them in the formalism. For example, for the 
simple case of only one extra neutral heavy lepton, the effects on a neutrino electron 
scattering off electrons will be a global factor in the Lagrangian, due to the non-unitarity 
of the mixing matrix. In this case the corresponding NSI parameters for an electron 
neutrino experiment will be given as 

£ e ee = -9 l sin 2 6>14, £ e e e = ~9 r sin 2 014, (9) 

where 6 14 is the mixing angle between the light neutrino and the extra heavy fermion 
and gL,R are the SM coupling constants for the neutrino electron scattering process. 

2.3. Additional scalars 

Despite the NSI formalism preserves the V — A structure of the theory, it is also possible 
to consider the impact of scalar couplings. For example, if we consider the case of 
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low energy supersymmetry with broken R-parity [8?, ESI SSJ where one has trilinear L 
violating couplings of the form 


^ijkLiL i E i 


ki 


Kjk L iQj D t 


( 10 ) 


with L and Q super-fields that contain the usual lepton and quark SU( 2) doublets, E c , 
and D c super-fields that contain the singlets, and i,j,k the generation indices. These 
couplings give rise, for example, to the following four-fermion effective Lagrangian for 
neutrino interactions with d-quark 

£ eff = -2 V2G f E £ d a ^LaY^Lf}d R Yd R: (11) 

a,/3 


where we have, among others, flavor-conserving and changing NSI, given, respectively, 
by 


dR 

b fin 


y l_ A t-il 2 alld = y . 

i ^G F m\ L ^ Z? 4 j2G F m% L 


( 12 ) 


Here, mq jL denotes the masses of the squarks, while j = 1,2,3 stands for d R , s R , b R , 
respectively. This is just one example, but other NSI couplings can be studied in this 
context [ID]; moreover, constraints on generic scalar NSI can also be studied by using 
Fierz transformations i!H. 


3. NSI phenomenology in propagation 

In this section, we discuss the phenomenological impact of the presence of NSI in 
propagation. The effect of NSI can be present through the modification of the matter 
potential that can exist not only in the diagonal but also in the off-diagonal elements in 
the effective Hamiltonian. Before starting our discussion on NSI we will briefly review 
the current standard oscillation status. 


3.1. Standard neutrino oscillations picture 

Unless otherwise stated, the standard three-flavor picture of neutrinos is assumed, u e , 
and v T and corresponding anti-particles. In vacuum, the mixing of neutrinos is supposed 
to be described by the usual flavor mixing without NSI, 

3 

K) = ^2 U ai Wi) (<X = e,fJL,T), (13) 

2—1 

where U is the 3x3 matrix which describes the flavor mixing [92] of neutrinos. In this 
review, we use the standard parameterization found, e.g, in [8], 
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C 12 C 13 S 12 C 13 s i 3 e~ lScp 

-S12C23 - C 12 S 2 3Si 3 e lScp C12C23 - S 12 S 23 S 13 e lScp S23C13 
S12S23 - c l 2 c 23 si 3 e l5cp -C12S23 - s 12 c 23 s 13 e z5cp c 23 ci 3 


(14) 


where Sij = sin %, c l3 = cos 9 t] , and <5 cp is the Kobayashi-Maskawa [93] type CP phase 
for neutrinos. 

In addition to the mixing angles and CP phase, the mass squared differences of 
neutrinos, Am| = rnf — rrij with m; (z=l-3) being the neutrino masses, are the relevant 
parameters to describe neutrino oscillation. So far, all of these parameters have been 
measured with reasonably good accuracies except for the value of the CP phase and 
sign of Amjj (Am| 2 ). The positive (negative) sign of corresponds to the normal 

(inverted) mass ordering, often referred to as normal (inverted) mass hierarchy. 

Different groups have carefully studied the neutrino data and obtained accurate 
values for most of the three-flavor neutrino oscillation parameters j951 [5U . 157] . Their 
most important results are summarized in Table [lj where a reasonable agreement can 
be seen. 


Table 1 . Summary of the standard three-flavor picture parameters, as reported from 
three different groups, denoted as C (Capozzi el al [2SJ, second and third column, F 
(Forero et al |5B3, fourth and fifth column), and G (Gonzalez-Garcia et al [57], the 
last two columns). The parameter Am]]; has a slightly different definition in each 
case, being A= to§ — (m\ + m\)/2 for Ref. [95] > Am]], = m ‘l — mf for Ref. [96] 
and Am| ; = — m\ for normal hierarchy (NH) and A= m§ — m| for inverted 

hierarchy (IH) for Ref. |97J 


Parameter 

C [95] 

Best fit 3 cr range 

F [96] 

Best fit 3 cr range 

G [97] 

Best fit 3 a range 

Am^/10 - 5 eV 2 

7.54 

6.99-8.18 

7.60 

7.11-8.18 

7.50 

7.02-8.09 

A?n 2 ; /10 - 3 eV 2 (NH) 

2.43 

2.23-2.61 

2.48 

2.30-2.65 

2.457 

2.317-2.607 

-Am|/10 - 3 eV 2 (IH) 

2.38 

2.19-2.56 

2.38 

2.20-2.54 

2.449 

2.307-2.590 

sin 2 9 12 / 10 _1 

3.08 

2.59-3.59 

3.23 

2.78-3.75 

3.04 

2.70-3.44 

sin 2 # 23 /10 -1 (NH) 

4.37 

3.74-6.26 

5.67 

3.93-6.43 

4.52 

3.82-6.43 

sin 2 023/1O- 1 (IH) 

4.55 

3.80-6.41 

5.73 

4.03-6.40 

5.79 

3.89-6.44 

sin 2 d 13 /10 - 2 (NH) 

2.34 

1.76-2.95 

2.26 

1.90-2.62 

2.18 

1.86-2.50 

sin 2 6 * 13 /10 —2 (IH) 

2.40 

1.78-2.98 

2.29 

1.93-2.65 

2.19 

1.88-2.51 

5/° (NH) 

250 

0-360 

254 

0-360 

306 

0-360 

8/° (IH) 

236 

0-360 

266 

0-360 

254 

0-360 
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3.2. Neutrino evolution with NSI 


Phenomenologically, the evolution equation of neutrinos in the flavor basis in the 
presence of propagation NSI in unpolarized matter can be generically written as, 

?f 
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,( 16 ) 


where v cx = (u a \u(r)) ( a = e,/i,r) denotes the probability amplitude to find neutrino 
as u a at the position r, A ;J - = Am|/( 2 P), if being the neutrino energy. V/ = V/(r) = 
y/2 GfN f{r) where A/ (/ = e, w or d) denotes the fermion number density along the 
neutrino trajectory in matter. Note that V e (r) is the standard matter potential |16j 
which induces the usual MSW effect [93, 16] , 

Since NSI effects in propagation enter only through the vector couplings, must 
be interpreted as e^p = + e^. For simplicity, throughout this review, we consider 

the case where only d-quark has the propagation NSI with neutrinos, and write NSI 
parameters simply as e a p by omitting the fermion superscript. Note that the case of 
u-quark NSI is very similar in most cases to be discussed in this section because in the 
usual matter, N u ~ N d ~ 3 N e . 

In this review, unless otherwise stated, for definiteness, we use the following values 
of the mixing parameters as our reference values; sin 2 6 12 = 0.31, sin 2 #13 = 0.023, 
sin 2 6*23 = 0.5, A= 7.5 x 10 -5 eV 2 and (Am 2 ^ = 2.4 x 10 -3 eV 2 , and Sep = 0, which 
are consistent at 2 cr with the results obtained by the recent global analysis [951 |96l [97]. 

Eq. (15) defines the framework for neutrino propagation in matter with NSI. The 


parameters e a p (a, f3 = e,p,r ) describe the magnitude of NSI. The diagonal NSI 
parameters, e aQ (a = e,/i, r), could play a role similar to the terms of the standard 
MSW matter potential, or could be interpreted as the NSI induced mass squared 
difference, mimicking the ones that contain Am 2 , which could induce new resonance 
even if neutrinos were massless [181120] . On the other hand, off-diagonal NSI parameters, 
e a fs{a 7 ^ /3) could play a role similar to the mixing angle. Even if there is no mixing in 
vacuum, the flavor transitions u a —> vp can occur in matter due to the presence of the 
off-diagonal NSI [EEH1 CEE 2D] • The complex phases of the off-diagonal elements e a p could 
be a new source of CP violation, see e.g., mm- 

Currently, almost all the neutrino data are consistent with the standard three flavor 
scheme of massive and mixed neutrinos. Therefore, NSI, if they exist, is expected to 
manifest only as a subdominant effect. NSI in propagation has been constrained mainly 
by the oscillation data of solar and atmospheric neutrinos as well as neutrinos produced 
by accelerators. Reactor neutrinos do not constrain the propagation NSI because the 
matter effect is expected to be very small though they can constrain the detection NSI. 

Roughly speaking, for a given neutrino energy, and the matter density, p, the impact 
of propagation NSI in neutrino oscillation (modification of the standard oscillation due to 
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NSI) is essentially determined by the magnitude of the following dimensionless quantity, 

"2.4 x 10 ~ 3 eV 2 ' 


£-a/3 Vf 

^ 


0.1 x E a p 


E 


GeV 






P 


_g/cm 3 _ 


(16) 


where the NSI with d or u quarks are assumed; Am?- is the relevant mass squared 
difference in the corresponding oscillation channel, and the baseline L is assumed to 
be large enough, or A t] L > 0(1). Larger the value of r} a p, larger the impact of NSI in 
propagation. 


3.3. NSI for atmospheric neutrinos 


In this subsection we discuss the NSI effect for atmospheric neutrinos. The impact 
of NSI on atmospheric neutrinos have been considered by many authors, see e.g., 

Let us first consider the impact of NSI on the u /t — u T sector and assume that all 
the NSI parameters coupling to electron flavor neutrino, namely, e e /3 are zero. In this 
scenario, £ /iT and e TT — e /t/i can be constrained mainly by the higher energy samples of 
the atmospheric neutrino data as will be seen below. 

In the limit of A m^L/E —> 0, with the constant matter density approximation, 
and ignoring 9 i3 , the survival probability is expressed as (32, 42] 


p { v v GO 


1 — Piyfj, —> v T ) ~ 1 — sin 2 29 e g sin 2 


£ 


A31 L 
2 


(17) 


where 


sin 2 29 e s 


sin 20 23 ± 2?^ r | 2 

e 


(18) 


£ = yj | sin 26*23 ± ^Vrr\ 2 + {cos26» 23 + ~ Vtt)} 2 , (19) 

and the +(—) sign in front of r]^ T corresponds to the normal (inverted) mass ordering. 
For anti-neutrino channel, the sign of r] a p must be changed. 

For the atmospheric neutrinos, the sensitivity to the NSI parameters can be 
estimated by studying the muon neutrino and anti-neutrino survival probabilities for 
different zenith angles, cos 6 k. Fig. [l] shows the muon neutrino survival probabilities 
(for the normal mass ordering) for cos 6 L = —0.3 (left panels), —0.6 (middle panels) and 
—1 (right panels) for the cases without NSI (by solid lines) and with NSI (by non-solid 
lines). For this calculation, the neutrino evolution equation (15) was solved numerically 
(without ignoring neither Am^ nor $ 13 ) using the Earth matter density profile predicted 
in the Preliminary Reference Earth Model (PREM) model |98lj . 

Note that, by comparing the upper and lower panels, the dependence on the sign 
of NSI parameters for neutrinos and anti-neutrinos are opposite. Note also that, with 
a good approximation, the survival probabilities are invariant under the simultaneous 
transformation A 31 —>• — A 31 and e^ T —> In the limit of the 2 flavor approximation, 

what is relevant is only the relative sign of these quantities. As can be seen from Fig. 
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Figure 1. Survival probabilities of —>• is fl (upper panels) and —>• (lower 

panels) as a function of the neutrino energy for different zenith angle of incoming 
neutrinos, cos6 z = —0.3 (left panels), —0.6 (middle panels) and —1 (right panels). 
The corresponding distances traveled by neutrinos are indicated in the plots. The 
normal mass ordering was assumed. 


[lj the impact of NSI is small for lower energies, for E < 5 GeV, even for the case where 
neutrino pass through the center of the Earth. On the other hand, the impact of NSI 
for energy > 10 GeV, could be quite large for the NSI parameters considered in Fig. [l] 
and it is expected that these values could be disfavored or excluded. 

Here we quote the bounds on these NSI parameters obtained by the Super- 
Kamiokande collaboration [99], | e flT < 1.1 xlCT 2 and —4.9xlCT 2 < e TT — < 4.9 xlCT 2 
at 90% CL. More recently, by using the IceCube-79 and DeepCore data, authors of [1U01J 
obtained somewhat better bounds, |e MT | < 6 x 10 -3 and \e TT — e w | < 3 x 1CT 2 at 90% 
CL. 


For the v e — v T sector, besides the probability for u fl —> z/^, it is also useful to 
consider the —y v e case. The computation of these probabilities, shown in Fig. [2J was 
done in an analogous way to the case shown in Fig. [l] For this computation different 
NSI parameters have been considered: £ ee , e rr and e eT . It is important to notice that, in 
this case, the e eT parameters, could play a role similar to # 13 . Therefore, there is some 
impact on the —» v e channel as it is possible to see in the lower panels of Fig. [2j 

When £ ee , e tt and £ er are assumed to be simultaneously nonzero, it is known m 
that the allowed combinations of the NSI parameters are approximately given by the 
parabolic relation, 

3|e eT | 2 

£tt ~ - 

1 + 3£ e , 

This feature can be also confirmed in Figs. 8 and 9 in [99] which show the constraints 
on these NSI parameters obtained by the Super-Kamiokande atmospheric neutrino data. 
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cos0 =-0.3 (L = 3822 km) 


cos0 = -0.6 (L = 7645 km) 


cos0 =-l (L= 12742 km) 
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Neutrino Energy [GeV] 
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Figure 2. Survival probabilities of ly —► zy (upper panels) and zy —> v e (lower 
panels) as a function of the neutrino energy for different zenith angle of incoming 
neutrinos, cos0 z = —0.3 (left panels), —0.6 (middle panels) and —1 (right panels). 
The normal mass ordering was assumed. 

It is also pointed out in H that atmospheric neutrino data alone can not essentially 
constrain £ ee parameter. Therefore, in general one can obtain the allowed regions of e eT 
and e TT for given values of e ee as done in [43:, [99]. For example, from Fig. 9 of [99]. 
we see that for sin 2 #23 = 0.5, for £ ee = —0.5, 0, and 0.5, £ eT < , 0.08, 0.11 and 0.18, 
respectively, at 90% CL. 

3.4- NSI for accelerator neutrinos 

So far the bounds on propagation NSI from accelerator neutrinos mainly come from the 
zy —> zy and zy —)■ D M channels. For these channels, at first approximation, the relevant 
NSI parameters are £ MT , and £ TT . The zy —> zy and zy —> survival probabilities 
as a function of neutrino energy for the MINOS baseline, L = 730 km, are shown in 
Fig. i The computations were done without the presence of NSI (solid lines) and with 
£^ r = ±0.1 (dotted and dashed lines). As it is possible to see from Fig. [3] the impact 
of NSI for v and v channels are opposite. 

The bounds obtained by the MINOS collaboration jl lQl] . translated to the notation 
used in this review, can be stated as —0.067 < £ iir < 0.023 at 90% CL. 

3.5. NSI for solar neutrinos 

The most updated analysis of solar neutrinos in the context of the propagation NSI 
comes from m ■ For solar neutrinos, under the so called one mass scale dominance 
approximation, jAm^l —* 00 , the neutrino evolution can be effectively reduced to that 
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normal mass hierarchy (L = 730 km) sinD, 3 = 0.5 sirf20 13 = 0.089 S cp = 0 



Figure 3. v fl — > v tJ and —> ly, survival probabilities as a function of neutrino 
energy for the MINOS baseline, L = 730 km without the presence of NSI with and 
with NSI, e^t = ±0.1. The matter density was assumed to be constant, p = 3.2 g/cm 3 . 


of the 2 flavor system |102j . as follows, 
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with u' T being decoupled from the system [102], and diagonal £p> and off-diagonal e^ 
NSI parameters are related to E a p as [73], 

£ D — C13S13R-G [e l ' 5cp (s'23 c e/t + C23 £ eT )j — (1 + Si 3 )c 23 S 23 R- e ( e /iT) 
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S 23 


C 23 £ *llt + c 23S23(^rr — £ 


fifi 


In this approximation, the survival probability is given as 

P{y s -> v s ) — Sj 3 + 1 yj. 


(23) 
)]• (24) 

(25) 


where /h,,!c( ,y( j i s calculated for the effective 2 flavor system described by (21). 

According to U3), the bounds on these parameters are —0.25 < Ep, < —0.02 and 
—0.14 < e n < 0.12 at 90% CL assuming NSI with d-quark. 

In Table [2] we show the summary of the bounds on the propagation NSI. 


4. NSI phenomenology in detection 

Several experiments have been devoted specifically to measure with precision the 
neutrino interaction with quarks and leptons. They are performed at very short 
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Table 2. Constraints on the matter (propagation) NSI parameters at 90% C L. for 
the interaction of neutrinos with d type quark. The other NSI parameters are set to 
zero. 


NSI parameters 

Bounds 

Reference 

~m _ _m 

fc ee 

(0.02, 0.51) 

m 

_m _ _m 

^TT 

(-0.01,0.03) 

m 

~m _ ~m 

‘-'TT ULl 

(-0.049,0.049) 

[eg 

£.m _ 

b TT 

(-0.036,0.031) 

m 

~m 

(-0.09,0.04) 

m 

~m 

^ [IT 

(-0.01, 0.01) 

m 

~m 

b [lT 

(-0.011,0.011) 

pg 

~m 

b [lT 

(-6.1 x 10~ 3 , 5.6 x 10“ 3 ) 

[100) 

t-'er 

(-0.13,0.14) 

m 

C (for oS = 0.50) 

(-0.05,0.05) 

p9] 

C (for 0.50) 

(-0.19,0.13) 

[99] 


baselines, avoiding effects coming from the standard oscillation. These measurements 
allow us to test the validity of the interactions described by the Standard Model, and, 
therefore, they could be a basis to search for new physics beyond SM. 

For non-oscillation experiments, NSI can be constrained by comparing the measured 
cross sections with that predicted by SM for the interaction of the neutrinos with the 
corresponding target. Most of these experiments record fewer events than oscillation 
experiments. On the other hand, they are independent of the mixing parameters; 
therefore, the cross section measurements, in general, do not suffer from the uncertainties 
of the oscillation parameters. Moreover, non-oscillation experiments are sensitive to 
axial couplings, a coupling that is absent in propagation NSI effects. 

Here we will review different experiments that constrain NSI through detection. 
We will start by considering the neutrino interactions with electrons and, afterwards, 
we will review its interactions with quarks. 

We will illustrate the phenomenology involved in this type of experiments by 
considering the specific case of the electron anti-neutrino scattering off electrons. In this 
case, the differential cross section, including the corrections coming from the Lagrangian 
shown in Eq. (JT|) , will be given by 


da 

d T e 


2G 2 , 


,m P 


7T 


[(dR + £ 


R \2 


+ ^2 \ £ ae\~ + j (9L + £ee) 2 

a^e \ 
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( 9 L + £ee) (dR + £ f e ) + V. 


,n 


m f 


Ijl 

El 


• (26) 


Here, m e is the electron mass, T e = E e — m e (with E e being the total electron energy) 
stands for the electron recoil energy, and E v is the anti-neutrino energy. The Standard 
Model couplings, at tree level, are defined as ^ = 1/2 + sin 2 6\y and gn = sin 2 # w . 
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We can see from this expression that flavor changing NSI parameters (s^ e R and £ R,R ) 
will only give quadratic corrections while the flavor diagonal ones could give linear 
corrections. 



T e (MeV) 

Figure 4. Averaged differential cross section for the electron anti-neutrino scattering 
off electrons for the SM case (black solid line), for a flavor changing NSI (blue dashed 
line), and for a flavor conserving NSI (green dashed dotted line). The reactor anti¬ 
neutrino flux has been considered in order to integrate the anti-neutrino cross section 
over the appropriate neutrino energy range. 

This is illustrated in Fig. [4j where we show the differential cross section for anti¬ 
neutrino electron scattering, u e e~ —> u e e ~, averaged over a typical anti-neutrino reactor 
spectrum [10311104] . The plot is given in terms of the electron recoil energy, T e , for an 
energy window relevant for an anti-neutrino detector such as TEXONO [1051 106J. In 
the plot, the prediction for the SM cross section is shown, as well as that for the NSI 
one. For both flavor changing and flavor conserving NSI, the same negative value of 
the parameters are used; this illustrates how flavor diagonal NSIs have more impact on 
detection signals than flavor changing parameters. 

The NSI parameters for this reaction can be constrained by considering, for 
example, the data from the TEXONO collaboration, which use P e e scattering as 
the detection signal. We have updated the analysis reported by the TEXONO 
collaboration [105j . including the new predicted spectrum [ 103 . 1104j and radiative 
corrections [107j in order to obtain new constraints for these parameters. We have 
also combined the results of this analysis with the constraints coming from the v e e 
scattering measurements reported by the LSND collaboration [1'08] . By combining these 
two experiments we can obtain stronger bounds both on left and right NSI parameters, 
taking advantage of the different chirality of both neutrino experiments. The result of 
this new analysis is shown in Fig. [ 5 ] for the diagonal parameters e^. r , and is also shown 
in Table [3] along with other current constraints. 

Previous works involving the neutrino scattering off electrons obtained constraints 
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£ 

ee 

Figure 5. Allowed region, at 90 % CL, for diagonal NSI parameters, e^ R , from 
a combined analysis of TEXONO reactor anti-neutrino and LSND neutrino electron 
scattering off electrons. 

from a combined analysis of reactor neutrino experiments. There are different 
experiments that involve an electron (anti)neutrino flux, such as reactor neutrinos 
(TEXONO [105]. MUNU jl()9j . Rovno [110j . Krasnoyarsk [lTlj . Irvine [ TI2] ) and 
accelerator neutrinos (LSND [108] and LAMPF |113j ). 

We show in Table [3] the summary of the constraints from these analysis. We prefer 
to show in this table results obtained by different research groups since the analysis may 
have different assumptions that can be important in the interpretation of the parameters. 
These results consider either one or two parameters at a time, while all other parameters 
equal to zero; in some cases, especially for flavor diagonal couplings, the correlation 
between these two parameters is important. 

We also show in Table [3] the constraints obtained from muon (anti)neutrino fluxes, 
based on the results coming from the CHARMII experiments [ITS], Although there is 
no man made tau neutrino sources, it is possible to constrain these interactions if one 
considers the LEP measurements of the process [33] e + e _ —>- uhy where tau neutrinos 
appear as part of this inclusive reaction, or to consider the solar neutrino flux that also 
includes a tau neutrino component [60]; such constraints are also shown in Table [3] 

In order to get constraints on the NSI of neutrinos with d type quarks, it is 
necessary to study experiments such as CHARM, CDHS [ IT? . 1118] and, more recently, 
by NuTeV [ 119] . They have measured the cross section for the scattering of electron and 
muon neutrinos off quarks. For the case of NuTeV, there have been a long discussion 
about a discrepancy of the measured cross section with the SM prediction. After the 
revaluation of the predicted sea contributions from the c quark, it has been possible 
to solve this puzzle [1201 1121] ; currently NSI suggested by the NuTeV experiment are 
considered as consistent with zero. 

Table [djshows the summary of the constraints for the d quark NSI coming from these 
experiments. Constraints coming from charge lepton flavor conversion, such as /i —> eq 
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Table 3. Constraints on the detection NSI couplings at 90% C L. for the interaction 
of neutrinos with electrons 



one parameter 

two parameter 

~eL 

fc ee 

c-eR 

‘-'ee 

(-0.021,0.052) [60] 
(-0.07,0.08) [114J 

(-0.08,0.09) [T15] 

(-0.02,0.09) [68] 
(-0.11,0.05) [68] 

(-0.036,0.063) [60] 
(-0.10,0.09) [115] 

peL 

c-eR 

(-0.03,0.03) m 
(-0.03,0.03) jlH! 

(-0.03,0.03) [Mj 
(-0.03,0.03) pS] 

(-0.033,0.055) [5^ 
(-0.040,0.053) [M] 

~eL 

c rr 

p eR 

c Tr 

(-0.16,0.11) [60] 

(-0.46,0.24) [53] 
(-0.25,0.43) [53] 

(-0.51,0.34) [68] 
(-0.35,0.50) [68] 

(-0.16,0.11) [63 
(-0.4, 0.6) pS] 

c.eL 

fc e/i 

peR 

^efi 

(-0.19,0.19) [1T4J 

(-0.13,0.13) pS] 
(-0.13,0.13) [53] 

(-0.53,0.53) [33] 

(-0.53,0.53) pS] 

p eL 

^er 

c-eR 

t-'er 

(-0.4, 0.4) [33 
(-0.28,-0.05) and 
(-0.19,0.19) [114J 

(-0.33,0.33) [M] 

( 0.05, 0.28) [53] 

(-0.53,0.53) PS] 

(-0.53,0.53) [33] 

P-eL 

V 7- 

F eR 

^ [IT 

(-o.i, o.i) m 
(-o.i,o.i) m 

(-o.i, o.i) m 
(-o.i,o.i) m 

(-0.53,0.53) [33] 

(-0.53,0.53) [33] 


or fi —> e, have not been considered here. They always involve, at some level, a one loop 
dressing of the neutrino vertex; therefore, they will always be model dependent. The 
readers interested in such constraints can see, for example, those reported in |122l 165] . 

Table 4. Constraints on the detection NSI couplings at 90% C L. for the interaction 
of neutrinos with quarks 


NSI parameters 

Bounds Ref. 

p dL 
^ee 
pdR 
fc ee 

(-0.3, 0.3) [33 
(-0.6, 0.5) [33 

pdL 

pdR 

(-0.005,0.005) [S3 
(-0.042,0.025) [S3 

pdL 
^ fie 
pdR 
fie 

(-0.023,0.023) [S3 
(-0.036,0.036) [S3 

pdL 

b er 

pdR 

t-’er 

(-0.5, 0.5) [33 
(-0.5, 0.5) [33 

pdL 

pdR 

(-0.023,0.023) [S3 
(-0.036,0.036) [S3 


5. Future Prospects 

There are several experimental proposals that plan to improve the current knowledge of 
neutrino properties. Therefore, there is plenty of room to improve the sensitivity to NSI 
in the near future. Some of these experimental set-ups are discussed below, showing the 
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future perspectives for different types of experiments. Again, the discussion is divided 
into propagation and detection NSI. 

It is important to notice that, besides the need for more restrictive bounds on NSI 
parameters, it is also necessary to solve the possible confusion between standard and 
non-standard parameters. As has been stated in the past [36], it is possible to have 
a confusion in neutrino oscillation experiments between NSI parameters and standard 
mixing angles, especially 0 13 ; recently, this subject has been discussed in the context of 
solar [123 ] and reactor neutrinos [78]. The significant progress in improving the precision 
on 8 13 will strongly restrict this possibility in the near future. 

Another important topic in this direction is that of the robustness of the solar 
neutrino data against NSI. It might be possible that large NSI effect give rise to a dark- 
LMA solution without contradicting any current experimental result H23. This solution 
has persisted as a plausible picture [El EH] • A recent study on the future combined data 
of JUNO [ 125] . RENO-50 [ 126] and NOvA [ 127] has discussed the perspectives to exclude 
this solution [128]. 

5.1. perspectives for NSI in propagation 

Different experimental set-ups (proposed to increase the precision for the determination 
of the neutrino oscillations parameters) have been considered in the last years. The 
need for a better knowledge of the Kobayashi-Maskawa type CP phase in the lepton 
sector as well as the mass ordering is certainly a major motivation for these proposals. 
Currently ongoing experiments such as T2K and NOvA may improve somewhat the 
current bounds on some NSI parameters but probably not so much, especially for that 
coming from the v e (D e ) appearance mode due to relatively small statistics. 

Hyper-Kamiokande [HE OSD] is an interesting proposal that expects to improve 
the sensitivity to NSI, by using the atmospheric neutrino data [ 1311 132 . 133] ; their 
expectations for the normal hierarchy case are particularly appealing. On the other 
hand, for the case of the LBNE [134] and LBNO [135] proposals, the expected sensitivity 
to NSI is also encouraging, especially for the flavor changing case of e /ir and e efl 1136 . 

Finally, important constraints are expected from the IceCube Deep Core and 
PINGU experiments. These are extensions of the IceCube experiment focused on a 
lower energy range. In this case, the expectations to constrain the flavor changing NSI 
parameter e^ T could reach the one percent level. For the flavor diagonal case, it might 
be possible to obtain information about the elusive parameter e TT [100, 13711138] . 

5.2. perspectives for NSI in detection 

The future neutrino oscillation experiments will also be sensitive to the NSI parameters 
through a detection effect. For instance, for the case of the proposed JUNO [ 125] and 
RENO-50 [ 126] an improvement to the constraints on e e/i and e er is expected [77] . 

For the case of the interaction of electron neutrinos with electrons, both 
ISODAR [139] and LENA 11 Oj proposals could give complementary information if 
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both proposals are done in the future. In the case of ISODAR, it is proposed to 
use an intense anti-neutrino 8 Li source, with an anti-neutrino energy ranging up to 14 
MeV, in combination with the KamLAND liquid scintillator [139]. The LENA proposal 
plans to use a neutrino Chromium source, providing a monochromatic neutrino flux of 
energy, E u = 0.747 MeV, located at the top of a 100 kTon liquid scintillator cylindrical 
detector m- These are not the only proposals for neutrino electron scattering, but 
they illustrate the future potential of this experiments. 

The use of either a neutrino or anti-neutrino source leads to a better determination 
of the left or right-handed couplings, respectively. Therefore, if both neutrino and 
anti-neutrino experiments are done in the future, there could be a good room for the 
improvement of the NSI parameters. 

We illustrate this by showing, in Fig. [6j the expected sensitivity for the case of 
a neutrino artificial source in combination with the proposed LENA detector as has 
already been calculated in [141]. where an expected total number of 1.9 x 10 5 neutrino 
events and a 5 % systematic error was considered. The case of an anti-neutrino source 
is also shown in Fig. [6] In this last case the analysis developed in [13~9] has been closely 
followed. The expected result from the combined analysis of both future experiments 
is shown by the region filled by the magenta color. We show in the same figure one 
of the current constraints on NSI, coming from the solar neutrino analysis [BO]. It is 
possible to see that there is room for improving these constraints by almost one order 
of magnitude, especially if both experiments are realized. 



R 

8 

ee 


Figure 6. Expected sensitivity to NSI parameters e^ ,R from ISODAR and LENA 
proposals. The result of a combined analysis is also shown. 

Regarding the NSI of neutrinos with quarks, there are severals proposals that 
could improve current bounds. An interesting proposal is that of neutrino coherent 
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scattering off nuclei |142j . After seminal works on the construction of these type of 
detectors [1431 1144] , there was a renewed interest in the previous decade |l~4'5l. 446] . At 
the same time, the sensitivity to NSI parameters and new physics searches was also 
noted [ 147 ; 1491] . Different experimental set-ups have already been considered, either 
using a reactor anti-neutrino flux [145] . and spallation source [1481 449 . H5D] , beta 
beams (EH [152], or pion decay [ 153 ; 11541 155 ]. These experiments could have an 
excellent sensitivity to the NSI. In particular, the TEXONO proposal has been studied 
in the past, using either a 76 Ge or 28 Si as a detector; an improvement of even one order 
of magnitude could be achieved in this case m • Other possible nuclei have also been 
studied [122] such as 48 Ti and 2 'Ah 

6. Conclusions 

Neutrino experiments have shown the existence of a new sector beyond the Standard 
Model because of the experimental evidences of nonzero neutrino masses that is already 
part of the current knowledge on particle physics. The mass and mixing of lepton 
sector turned out to be non-trivial, very different from that of the quark sector. 
While our knowledge on neutrino properties are continuously improving, the theoretical 
explanation of the neutrino mixing and the neutrino mass pattern is still an open 
question. 

In this context, possible presence of the non-standard interaction of neutrinos and 
its impact was discussed in this brief review from a phenomenological point of view, 
describing the current status on the search for new physics coming from NSI. 

So far there is no experimental evidence or indication of the presence of NSI and 
there exist only the constraints, which are summarized in Tables [2j [3] and [4] where we 
show the lower and upper bounds of the NSI parameters. Ongoing as well as proposed 
near future neutrino experiments are expected to improve considerably the NSI bounds 
or may indicate the presence of NSI. 
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